ApJ submitted 20 March 2011, accepted 10 July 2011 

Preprint typeset using 1^1^^ style eniulatoapj v. 11/26/04 



O 
(N 



o 
u 

6 



> 

O 



% 



QUASISTEADY CONFIGURATIONS OF CONDUCTIVE INTRACLUSTER MEDIA 

G. M. VoiT^ 

ApJ submitted 20 March 2011, accepted 10 July 2011 

ABSTRACT 

The radial distributions of temperature, density, and gas entropy among cool-core clusters tend to 
be quite similar, suggesting that they have entered a quasi-steady state. If that state is regulated by a 
combination of thermal conduction and feedback from a central AGN, then the characteristics of those 
radial profiles ought to contain information about the spatial distribution of AGN heat input and the 
relative importance of thermal conduction. This paper addresses those topics by deriving steady-state 
solutions for clusters in which radiative cooling, electron thermal conduction, and thermal feedback 
fueled by accretion are all present, with the aim of interpreting the configurations of cool-core clusters 
in terms of steady-state models. It finds that the core configurations of many cool-core clusters 
have entropy levels just below those of conductively balanced solutions in which magnetic fields have 
suppressed electron thermal conduction to ~ 1/3 of the full Spitzer value, suggesting that AGN 
feedback is triggered when conduction can no longer compensate for radiative cooling. And even 
when feedback is necessary to heat the central ^ 30 kpc, conduction may still be the most important 
heating mechanism within a cluster's central ^ 100 kpc. 

Subject headings: cooling flows — galaxies: clusters: general — X-rays: galaxies: clusters 



1. INTRODUCTION 

The cores of galaxy clusters have lessons to teach 
about accretion of hot intergalactic gas onto massive 
galaxies and the feedback processes that limit star for- 
mation within those galaxies. Many of the most mas- 
sive galaxies in the universe can be found at the cen- 
ters of galaxy clusters. A cluster's central galaxy is 
often the brightest cluster galaxy (BCG), and the ma- 
jority of central galaxies exhibit little or no star for- 
mation. Minimal star formation is understandable for 
central galaxies in which the hot intracluster gas can- 
not radiate its thermal energy within a Hubble time but 
is harder to understand in clusters with shorter central 
cooling times. Star formation is generally seen in the cen- 
ters of gala xy clusters only when the central cooling time 
is < 1 Gvr (iMc Namara fc 0'Conncll"199a'CaIdiel et al.l 
ll998HRaiffertve t al. 2008; BildfeU et al. 2008). And even 
then, the star formation rate is usually < 10% of the rate 
one would naively infer fro m the apparent co oling rate 
of the intracluster medium (jO'Dea et al.l 120081 ). 

Feedback from a central active galactic nucleus (AGN) 
is the favored mechanism for suppressing cooling and 
star formation at the ce nters of these clusters (see 
IMcNamara fc Nulsenll2007l . and references therein) . Rel- 
ativistic plasma flowing from the AGN apparently evac- 
uates large cavities within these clusters and in at least 
some cases drives sh ocks into the surrounding X-ray 
emitt ing medium fe.g.^ Forman et al1l2007t [Randall et al.l 
1201 Ih . The energy introduced is comparable to that ra- 
diated from the ambient hot gas, making t he hypothesis 
of AG N feedback seem quite plausible (e.g. iBirzan et al.l 
[200l . 

Yet, many questions about the AGN feedback mech- 
anism remain unanswered. Two of the most pressing 
are: (1) how is the energy output of the AGN so well 
tuned to match the cooling rate of the surrounding hot 
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gas, and (2) how does the AGN manage to distribute 
its energy throughout the core so that heating nearly 
balances cooling everywhere? These questions arise be- 
cause the AGN accretion rate would seem to depend 
on conditions in a neighborhood only a few tens of 
parsecs in size, far smaller than the '-^ 100 kpc clus- 
ter core that must be heated. Furthermore, the en- 
tropy gradients of cluster cores with short central cool- 
ing times tend to increase gradually from 10 — 100 kpc. 
A large AGN outburst that dumped most of its en- 
ergy within a ~ 10 kpc region would overheat the cen- 
ter of the cluster, pro ducing an entropy inversion (e.g. 
IVoit fc Donahuel[2005D . Instead, the energy source that 
offsets coohng somehow balances the cooling rate out to 
~ 100 kpc without disrupting the gradually increasing 
radial entropy gradient. Some implementations of AGN 
feedb ack in galaxy-cluster simulations appear promisin; 
(e.g.. iBriiggen fc Sca nnapiec ol 1200^ iDubois et al.l [201 
iFabian et al.ll2010l :lMcCarth v et al.ll201lD . but the prob 
lem is far from being solved. 

A closer look at the data suggests that clusters with 
central cooling times < 1 Gyr do indeed achieve a quasi- 
steady balance between heating and cooling that some- 
times allows star formation to proceed at moderate rates 
of 1 — lOMoyr"-'^. Gas temperatures at ~ 10 kpc radii 
in clusters with short central cooling times are 2-3 times 
cooler than those at ^ 100 kpc, which has led such clus- 
ters to be called "cool-core" clusters. Gas densities de- 
cline by a factor ^ 10 over the same range in radius. And 
even though gas density and temperature proflles in these 
systems differ in normalization, their proflles of entropy 

{K — Pp^^/^ (X kTue ) as a funct i on of radius are 
all quite similar ("Donahue et al.l [20061 ICavagnolo et al.l 
2009; Sanderson et al. 2009). A cluster's entropy dis- 
tribution is of interest because it determines the tem- 
perature and de nsity proflles of a cluster in hydrostatic 
equihbrium (see lVoit et al.l (|2002D and § 2). Mergers and 
accretion shocks driven by gravitational structure for- 



mation produce a family of clusters with self-similar en- 
tropy profiles, and deviations from those gravitationally 
generated profiles reflect the actio n of non-gravitat ional 
cooling and heating processes Ce.g. lVoit et al.|[2003l) . 

The purpose of this paper is to explore why one partic- 
ular set of density, temperature, and entropy profiles is 
favored among cool-core clusters. Its mode of exploration 
is to outline the properties of various families of steady- 
state solutions for hot gas within massive, spherically- 
symmetric dark-matter halos and compare them with 
observations of cool-core clusters. While the outer re- 
gions of clusters are certainly not in a steady state be- 
cause of ongoing accretion and merger events, the inner 
regions are likely to be closer to a steady state because 
the relevant timescales for sound waves, gravitational in- 
teractions, and radiative cooling are all ^ 10 Gyr at 
r < 100 kpc. 

Our exploration proceeds as follows. Section [2] sets the 
stage by presenting some useful expressions for the struc- 
ture of clusters in hydrostatic equilibrium and demon- 
strating that the temperature profiles of cool-core clus- 
ters are determined primarily by the shape of the under- 
lying gravitational potential, as long as the entropy of the 
central gas is sufficiently low. Section [3] then looks at the 
entropy profiles expected of cluster cores in which cooling 
is uncompensated by feedback. These pure-cooling pro- 
files turn out to be nearly power-law in form, K(r) ex r" 
with a ~ 1.2 at large radii. However, the implied cooling 
rates of clusters in these configurations are well known to 
be far larger than the star-formation rates seen in BCGs. 
Some form of heating must be compensating for radia- 
tive cooling, and we would like to know how the entropy 
profiles of clusters respond to that heating mechanism. 

Section 3] brings conduction into the mix. Con- 
duction in the context of cool-core clusters has a 
long and checkered history. It has been pro- 

posed several times as a mechanism with the po- 
tential to bala nce radiati ve cooling in cluster cores 
(iTucker fc Rosn cr 1983; B ertschinger fc MeiksinI 119861: 
iNaravan &: Medvedevii2001i) . However, a cluster in which 
electron thermal conduction b alances radiative cooling 
is in an unstable configu ration ( Bregman fc Davidlll988l : 
lSokeHl2003HGuo et al.ll2008.1 . Furthermore, some cluster 
cores have cooling rates too large to be balanced by elec- 
tron thermal conduction, even if it operates at its max- 
imuin rate, unfettered by magn etic fields (|Voigt et"al] 
I2002t iZakamska fc Naravanir2003D . The inability of con- 
duction to suppress cooling flows has led some aut hors to 
argue that it is unimportant in cluster cores (e.g.. lSok"eil 
120101) , and there are certainly sharp cold fronts in cluster 
cores across which conduction must be highly suppressed 
(jMarkevitch fc Vikhlinirj|2007l) , presumably by magnetic 
fields running parallel to the front. But magnetic fields 
cannot suppress conduction simultaneously in all direc- 
tions, and recent work has shown that MHD instabilities 
driven by anisotropic conduction may acti vely sculpt the 
magnetic field geometry i n cluster cores (|Parrish et al.l 
[2008: Parrish fc Quataerti[2008: McCourt et al.ll2010[ ). 

This paper does not investigate the fascinating astro- 
physics that can result from anisotropic conduction. In- 
stead it adopts a phenomenological approach in which 
conductive heat flow runs parallel to temperature gradi- 
ents, modulo a scalar suppression factor fc to account for 
the presence of magnetic fields. If turbulent velocities in 



the intracluster medium are sufficiently large, they will 
overpower buoyancy- driven MHD instabilities and ran- 
domize the magnetic field, leading to a mean local sup- 
pression f actor (fc) ~ 1/3 relative to a randomly chosen 
direction (Parr ish et all 120101 : iRuszkowsk i fc OhI |2010[ ) . 
We will therefore adopt fc = 1/3 as a fiducial value for 
the scalar suppression factor. 

Even though thermal conduction alone cannot stably 
balance radiative cooling, it remains interesting because 
it might be g overning the triggering of AGN feedback in 
cluster cores. iDonahue et al.l (|2005f ) raised this possibil- 
ity in an analysis of clusters with central cooling times 
much less than a Hubble time but with no obvious evi- 
dence for central star formation or AGN feedback. Those 
clusters had core entropy levels ~ 30 — 50 keVcm^, corre- 
sponding to cooling times ~ 1 Gyr, whi ch could in prin- 
ciple be offs et by thermal conduction. ICavagnolo et al.l 
()2008l 120091 ) confirmed this relationship between core en- 
tropy and feedback in a much larger sample: Clusters 
with core entropy < 30 keV cm^ often show evidence for 
AGN feedback while those with higher core entropy lev- 
els generally do not. There are some clusters and groups 
with powerful central radio sources that appear to have 
core entropy levels exceeding ^ 30keVcm^, but closer 
inspection reveals that these apparent exceptions have 
coronae of denser, lower-entropy gas i n a kilopar sec-scale 
region surrounding the central AGN ()Sunll200"9l) . 

Several studies have shown that thermal conduction, 
perhaps in concert with AGN feedback, can plausi- 
bly ex plain the observed bimodalit y in core behavior 
(IVoit et al. 2008; G u o fc Oh 2009; Sanderson et all l2009l: 
iParrish et all 120101: IRuszkowski fc OhI |2010D " Objects 
with core entropy < 30keVcm^ have central cooling 
rates that cannot be balanced by conduction, and AGN 
feedback is therefore necessary to stabilize cooling. In 
objects with greater core entropy, conduction is more 
efficient than cooling, as long as fc ^ 0.2, and a mod- 
est amount of intracluster turbulence may be necessary 
to ra ndomize the magnet ic fields and keep fc sufficiently 
large (jParrish et al.|[20Tot ) . Feedback from an AGN is not 
essential to counteract cooling in these objects because of 
the lower gas density in their cores, and both conductive 
and dynamical heating are capable of boosting the core 
entropy until the central cooling time exceeds several Gyr 
(jConrov fc Ostrikeij[2008l : IParrish et al.ll2010D . 

Section |4] addresses these issues by solving the steady- 
state equations for intracluster-medium configurations 
in which electron thermal conduction balances radiative 
cooling. While conductively balanced solutions for clus- 
ters can differ greatly in density and temperature, they 
are restricted to a narrow locus in the K-r plane and 
can be approximated with a simple power-law expres- 
sion. This locus is insensitive to cluster temperature for 
objects warmer than ~ 3 keV but depends somewhat on 
fc- Comparing the locus of conductive balance with the 
entropy profiles of actual clusters shows that the entropy 
slopes of cool-core clusters tend to break when they cross 
the conductively-balanced locus. At large radii, the K(r) 
profiles of all clusters lie above the conductively balanced 
locus and have slopes consistent with gravitational struc- 
ture formation. At small radii, the entropy proflles of 
many cool-core clusters lie at or slightly below the crit- 
ical locus for conductive balance, with slopes that track 
the conductively balanced conflguration. This behavior 



suggests that AGN feedback switches on when thermal 
conduction can no longer compensate for radiative cool- 
ing. 

Section [5] adds the element of thermal feedback. When 
conduction is present, steady-state solutions with inward 
accretion are not allowed unless cooling can outpace con- 
ductive heating. This section shows how steady-state so- 
lutions including AGN heating modify the conductively- 
balanced configuration, allowing lower entropy levels to 
persist in regions where AGN heating exceeds conduc- 
tive heating. It also compares steady-state models with 
cooling, conduction, and thermal feedback with obser- 
vations of real clusters. These comparisons show that 
conduction may be a more important heat source than 
AGN heating in the cores of many cool-core clusters and 
that AGN heating is usually necessary to counterbalance 
cooling within only the central '^ 30 kpc. 

Section |6] concludes by summarizing these results. 

2. HYDROSTATIC EQUILIBRIUM 

Everyone should know by now that the hot gas in a 
galaxy cluster is not in strict hydrostatic equilibrium. In- 
completely thermalized gas motions stimulated by merg- 
ers may provid e 5%-20% of the pressure support (e.g., 
iLau et aLll2009[ ). And heating or cooling of gas in the 
central region could drive either outflow or inflow. How- 
ever, because these motions are generally subsonic, hy- 
drostatic equilibrium is likely to be an adequate approxi- 
mation for the majority of relaxed-looking clusters. This 
section presents some useful formulae that describe hy- 
drostatic cluster structure and illustrate how the tem- 
perature structure of a galaxy cluster depends on the 
shape of its gravitational potential well and the entropy 
profile of the intracluster gas. In particular, these solu- 
tions demonstrate that clusters with cool gas at small 
radii do not necessarily contain gas that is actively cool- 
ing. The positive radial temperature gradients seen in 
cool-core clusters are primarily a reflection of the halo's 
gravitational potential shape and demonstrate that th e 
potential is close to the NFW form (jNavarro et al.llT997l ). 

Observational limits on the gas cooling rates in galaxy 
cluste rs have e stablishe d that they are not extremely 
large ([Peterson fc Fabianll2006[ ) . And as long as the cool- 
ing rate within a galaxy cluster is moderate, it will not 
drive the core away from hydrostatic equilibrium. For 
example, the steady-state inflow velocity at radius r is 
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where M is in MQyr~^, r^pc is r in kpc, p is the 
gas density, and rp is normalized to a value typical of 
the in ner 100 kpc of cool-core clusters (jVoit fc Donahue! 
l2005f) . Therefore, the inflow velocities induced by moder- 
ate rates of core cooling (< 10 Mq yr~^ ) should be highly 
subsonic in the regions resolvable with current X-ray tele- 
scopes. 

In the limit of negligible inflow velocity, the steady- 
state momentum equation reduces to one of hydrostatic 
equilibrium, which can be expressed as 

A p(7-l)/7 ^ _7^_GM_ 

dr 7 lO/ir'^ ^ ' 

for gas with equation of state P — Kp^ . Assuming that 
the gaseous contribution to Mr is small, one can directly 



integrate this equation for a given entropy profile K{r), 
yielding 
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where Pi, Ki, and Ti are the pressure, entropy, and tem- 
perature at some fiducial radius ri . Equation ([3]) is a gen- 
eralization of the pressure solution for p ower-law entropy 
profiles given bv ICavaliere et all (|2009( ). Analogous hy- 
drostatic solutions for the density and temperature pro- 
files follow directly from the relations p = {P/ KY^^ and 

These relationships can be used to clarify how a clus- 
ter's gravitational potential well and intracluster entropy 
gradient together determine the run of gas temperature 
with radius. Realistic configurations for gravitation- 
ally confined intracluster media have substantial pressure 
gradients for which P ^- at large radii. Equation ([3]) 
therefore implies 
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where r^, is a boundary radius at which P <C Pi and 

GMrprrip 1 2 

kT^ = Yj, = 2^"^P"c 

defines the characteristic temperature T^ and circular ve- 
locity Vc associated with the mass Mr within r. In other 
words, equation Q shows that the temperature of the in- 
tracluster medium at radius ri is determined almost ex- 
clusively by the depth of the potential well and the shape 
of the entropy profile exterior to that radius. It does 
not depend on the normalization of the entropy profile 
because multiplying the both the pressure and density 
profiles by the same factor C corresponds to multiplying 
the entropy profile by a factor C^~'~' without changing 
the temperature profile. 

Notice also that for a typical power-law entropy profile 
K (X r" with a ^ 1, the temperature at a given radius 
depends primarily on the value of T^ at slightly larger 
radii and on the value of a. The integral over T^ near the 
boundary radius carries less weight. For an isothermal 
potential, integrating over a power-law entropy profile 
gives Pi w [2(7 — l)/a]T0, demonstrating that shallower 
entropy profiles lead to greater gas temperatures within 
a fixed gravitational potential (see also Figure [Ij . 

The rising temperature gradients generally observed 
from ~ 10 kpc to ^ 100 kpc in cool-core clusters therefore 
refiect the fact that the mass profile within ~ 100 kpc 
rises more steeply than an isothermal Mr ex r profile. 
Consider the case of a gravitational potential consisting 
of an NFW halo with concentration ca relative to the 
radius ta, for which 
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■ CAr/rA 

In order to account for the gravitational infiuence at 
small radii of a large central galaxy, one can add an 
isothermal mass profile corresponding to the constant 
value T<^,BCG- Solid lines in Figure [1] show the form 
of T^ir) = T^,-NFw{r) + T^^bcg for ca = 4 and cen- 
tral galaxy mass ratio P0,bcg/T'0,nfw(^a) equal to 0.0 



and 0.3. All of these temperature profiles are scaled to 
Ta = Tcj,{rA), and in each case, T^ steadily rises out 
to radii ^ 0.3rA and becomes nearly constant at larger 
radii. 

Temperature profiles of intracluster media within these 
gravitational potentials tend to peak at radii smaller than 
the peak of T^, as one would anticipate from equation 
([4]). Dot-dashed lines in Figure [1] show T{r) for a = 1.2 
and the two potentials illustrated by the solid lines; long- 
dashed lines show the same thing for a = 1.0. The outer 
boundary condition in both cases is T{rA)/TA — 0.65. 
As with Tff,{r), a greater BCG contribution to the poten- 
tial leads to a greater temperature at small radii. Notice 
that the shallower power-law entropy slope of a = 1.0 
leads to a peak temperature ~ 20% greater than for 
a = 1.2 and that the presence of a BCG potential pre- 
vents the inner temperature from dropping much be- 
low ~ 0.3 Ta, even though the gas entropy approaches 
zero at small radii. Formally, the minimum value of the 
gas temperature is ^ T^^bcGj and will be greater if the 
central entropy does not go to zero, perhaps explaining 
why there is little spectral evidence for gas cooler than 
^ 0.3 Ta in cool-core clusters. 

Comparing these hydrostatic temperature profiles to 
those of real clusters shows that temperature profiles of 
cool-core clusters may be simply understood as resulting 
from intracluster gas with a nearly power-law entropy 
profile that resides in an NFW potential. Do tted lines in 
Figure [T] show fits by IVikhlinin et al.l ()2006l ) to four in- 
tracluster temperature profiles of clusters between 5 keV 
and 8 keV observed with Chandra. They are character- 
istic of other cool-core clusters and peak at radii and 
amplitudes similar to the peaks of the model profiles, in 
which the position of the gas temperature peak echoes 
the peak of T^ at larger radii. As long as K{r) oc r" 
with a ^ 1, the temperature structure simply responds 
to the shape of the potential well. 

Understanding the origin of these power-law entropy 
profiles is easy at large radii but less so at small radii. 
Numerical simulations show that gravitational struc- 
ture formation naturally produces ent ropy profiles wit h 
a ~ 1.1 — 1.2 outside olF cluster cores (jVoit et al.ll2005D . 
There is less consensus among simulations on what K{r) 
should be at small radii. When radiative cooling is not al- 
lowed, adaptive-mesh hydrodynamical methods typically 
prod uce more core ent r opy than smoot h-particle meth- 
ods (|Frenk et al.lll999t IVoit et al.ll2005| ). This discrep- 
ancy has been attribute d to the greater mix ing allowed 
in grid- based methods (jMitchell et al.l [20091) . However, 
cooling and feedback processes strongly influence the en- 
tropy profile in cluster cores and likely domi nate grav- 
itat ional heating at the smallest radii (e.g., IVoit et aLl 
120031) . Nevertheless, observations show that entropy pro- 
files of many cool -core clusters remain cl ose to if oc r^'^ 
at r < 100 kpc (|Cavagnolo et al.ll2009l ). The primary 
goal of the rest of this paper is therefore to understand 
what physical processes determine the slope of a clus- 
ter's entropy profile inside the core regions, because that 
is what determines the radial profiles of temperature and 
density. 

3. PURE COOLING 

Let us next consider quasisteady configurations in 
which pure radiative cooling determines the inner slope 



of the entropy profile. Moderate cooling will cause a slow 
inflow, and the steady solutions we derive in this section 
show that the entropy slope produced by cooling at small 
radii is coincidentally very similar to the K oc r^"^ slope 
produced by gravitational structure formation at larger 
radii. 

When cooling is present, we need to find solutions that 
satisfy both the equation of hydrostatic equilibrium and a 
steady-state entropy equation in which inward advection 
of heat balances radiative losses: 

^-^ = VA(T). (7) 

7 — 1 r d In r 

Here A(T) is a cooling function and we have adopted a 
sign convention in which negative v and M correspond 
to mass inflow. The physical processes determining the 
entropy gradient become a little clearer if we rewrite this 
equation as 

d\x\K _ r (7 - l)p^A _ tflow 

dlnr V P fcooi 

with 

P 

iflow = -r V icooi = 7 TT^T (9) 

(7-l)p-^A 

Approximate equality between the flow time and the 
cooling time at all radii therefore permits an approxi- 
mate power-law entropy-proflle solution for a particular 
M determined by the normalization of the entropy pro- 
file. 

Yet another rearrangement of the steady-state entropy 
equation from pure cooling helps to illustrate the rela- 
tionship between the entropy normalization and M: 
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In a pure-cooling configuration, one expects K -^ as 
r — J> 0. With this inner boundary condition and a given 
temperature profile T{r), one can integrate to obtain 
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(11) 



Notice that each pair of temperature and entropy profiles 
that solves equations ([3]) and pT|) can be extended into a 
family of solutions that share the same profile shape but 
whose entropy normalization is oc M'^^^'''/^. This rela- 
tionship arises because M is proportional to the square of 
the density scale in a pure-cooling solution, and entropy 
is proportional to p^~^ ■ Furthermore, if the temperature 
proflle is isothermal and 7 = 5/3, we obtain K (x r. 

Real clusters often have positive radial temperature 
gradients in their cores, making the effective power- 
law entropy slope of the pure-cooling solutions slightly 
steeper than linear. Figure [2] shows some steady pure- 
cooling solutions for the gravitational potentials in Fig- 
ure [1] and free-free cooling with A(T) oc T^/^. These so- 
lutions were obtained from equations ^ and (I10|) with 
the boundary conditions K/K^ = 1 and T/Ta — 0.5 at 
r/rA — 1, and M was adjusted so that K ^ a.s r ^ 0. 
They are therefore characteristic of any pure-cooling so- 
lution with realistic K and T at large radii and vanish- 
ing K at the center. They show that within the core of 
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Fig. 1. — Relationship between hydrostatic gas temperature and gravitational potential for power-law entropy profiles. Both panels show 
the same information, but the right-hand plot has a logarithmic radius scale to make the core structure more obvious. Two solid lines in 
each panel show the characteristic gravitational temperature T^ for halos consisting of an NFW profile with concentration ca = 4 plus an 
isothermal BCG contribution amounting to 71/),bcg/^</>,nfw(^a) = 0.0 and 0.3. Lines that are higher at small radii correspond to a greater 
BCG contribution. Two dot-dashed lines show the hydrostatic temperatures for gas in those halos, assuming K oc r^'^ and a boundary 
condition T jT/^ = 0.65 at r^- Long-dashed lines show the same thing for K oc r^ *-* . For comparison, dotted lines show temperature profiles 
(with r A = r5oo ) of the four clusters measured with Chandra by Vi khlinin et al.l II2005I . 12009 ) that have temperatures between 5 ke V and 
8 keV. This comparison shows that the temperature profiles of cool-core clusters may be simply understood as resulting from intracluster 
gas with a nearly power-law entropy profile residing in an NFW potential. 



a cluster, where cooling times can be short enough for 
quasi-steady cooling to be a possibility, the pure-cooling 
solution is similar to a power law with K oc r^'^. For 
the purpose of this plot, the entropy normalization fac- 
tor if A is arbitrary, because it is degenerate with the 
adjustable parameter M . 

The agreement in power-law slope between these pure- 
cooling solutions and the K oc r^-^ power law from grav- 
itational structure formation is purely coincidental but 
may help explain why the entropy profiles of cool-core 
clusters show no obvious break in slope where the intr- 
acluster cooling time equals the age of the universe. It 
is coincidental for the following reason: The requirement 
that iflow ~ icooi at all radii implies p oc r~'^/^(T/A)^/^ 
and K oc r(T^A)^/'^ for 7 — 5/3, meaning that the index 
of the resulting power-law relationship between K and r 
is determined in part by the temperature dependence of 
the cooling function, which has nothing to do with grav- 
itational heating. Because of this coincidence, there is 
no break in slope when a structure-formation solution at 
larger radii (> O.Ita) transitions to a quasisteady pure- 
coolin g solution at smaller radii (see lEttori fc Brigheniil 
I2OO8I for simulations demonstrating this behavior). 

The long-standing problem with such pure-cooling so- 
lutions for cluster cores is, of course, that they predict 
far too much cooling. Equation ()10p can be used to show 
how the excessive cooling rates of pure cooling solutions 



are related to the normalization of the entropy profile. 
Solving this equation for M with 7 = 5/3 and K oc r" 
yields 
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Intracluster plasma with a temperature exceeding 2 keV, 
has a cooling function approximately equal to 
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and the entropy profiles of clusters with cool cores 
tend to have K(r) ^ 1.6 x lO'^-^ergcm^ g~^/'^(r/lkpc) 
(e.g.. iCavagnolo et al.l 12009V In slightly less unfamiliar 
units, this representative entropy profile corresponds to 

kTue r~^ ~ 1.5keVcm^ kpc~^. Plugging this cooling 
function and entropy profile into equation (jl2[) gives the 
following relationship between the normalization of the 
entropy gradient and the mass fiow rate in the case of a 
steady pure-cooling solution: 
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Notice that the mass fiow rate in such a steady cooling 
solution is highly sensitive to the normalization of the 
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Fig. 2. — Entropy profiles for steady pure-cooling models. Solid 
lines show profiles for the two gravitational potentials from Figure 
1, and as in Figure 1, a greater BCG contribution leads to a higher 
line at small radii. Each entropy-profile solution has boundary 
conditions K/K^ = 1 and T/T/^ = 0.5 at r = r^, with M adjusted 
so that K —^ as r ^ 0. The dotted line shows a K oc r^'^ 
reference profile, illustrating the similarity of these steady pure- 
cooling solutions to that power-law profile. 



entropy profile. Reducing the entropy by 50% at all radii 
would increase the density and pressure by a factor of 
2"^/^ and the mass flow rate by a factor of 8. As discussed 
at the end of the next section, this strong dependence 
of steady-state mass inflow rate on entropy level should 
prevent core entropy profiles from declining too far below 
the critical profile at which feedback starts to operate. 

4. CONDUCTIVE EQUILIBRIUM 

Over the years there have been many efforts to find so- 
lutions in which thermal conduction alone balances cool- 
ing, bu t none has been complet e ly successful (see, for ex- 
ampl e , iTucker fc Rosned 119831; iBertschinger fc MeiksinI 
Il986t iNaravan fc Medvedevll20Qlt ). However, it is stih 
useful to understand the characteristics of such solutions, 
because conduction in combination with AGN heating 
can in principle lead to quasi-steady cluster configura - 
tions (jRuszkowski fc Begelma^ l2002t iGuo et all [20081 ). 
and the conductively balanced solution represents the 
limiting case in which AGN heating goes to zero. 

In the conductively balanced case, the steady-state en- 
tropy equation that must be solved in conjunction with 
hydrostatic equilibrium is 



r^ dr dr 



(15) 



where the thermal conduction coefficient k is a func- 
tion of temperature alone if electron diffusion dominates 
the heat flux. Typically, the density of the intracluster 
medium changes with radius much more rapidly than 
the temperature profile. It therefore makes sense to 
seek approximate solutions under the assumption that 
aT = dluT/dlnr is nearly constant with radius, in 



ar [1 + aria,^ + 1)] 
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(16) 



where «« = dluK/dlnr and the Field length (jFieldl 
119651 iBegelman fc McKedllQQOl ) is defined to be 



A, 




(17) 



Cool-core cl usters typica lly have qt '^ 0-3 at radii from 
10-100 kpc (|Voigt et al.l 2002; D onahue et al.H 2006') and 
if Spitzer conduction with a^ = 5/2 applies, then the 
left-hand side of equation [TBI is ~ 0.6. The run of den- 
sity with radius in a conductively balanced configuration 
is therefore determined by the condition that the Field 
length at a given radius be similar to the radius itself. 

The entropy profile corresponding to equation (J16l) 
takes on a particularly simple form if the Field length 
is determined by Spitzer conduction and free-free cool- 
ing. In that case. 



6 X 10~^T^/Vcergcm"is"^K"^/2 
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where fc is a scalar suppression factor depending on the 
magnetic field geometry. The Field length is then a func- 
tion of entropy alone: 



AF = 0.2kpc/, 
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ponahue et al.ll2005D .^ 

Figure [3] shows some representative solutions for con- 
ductively balanced intracluster media within the poten- 
tial well of a 5 X 10^* Mq cluster having a halo con- 
centration of CA = 4. These solutions were obtained 
by numerically integrating equations ([2]) and (|15|) with 
fc — 1/3 and they depend on the values of T and 
Uc at the inner boundary. One set of solutions (solid 
lines) has kT = 0.5kTA « 2keV and Uc = 0.45, 0.15, 
and 0.045 cm~'^ at the inner boundary. The other has 
kT = 0.3fcTA « 1.2 keV and the same set of inner den- 
sity values. For simplicity, equation ([13]) was used for 
the cooling function. 

These conductively balanced solutions have the follow- 
ing features. Greater inner density leads to a more el- 
evated temperature profile, because a steeper tempera- 
ture gradient is needed to balance the additional cool- 
ing. Greater inner temperature for a given inner density 
leads to a more elevated density profile, because a greater 
pressure is needed to confine the hotter medium. Among 
these solutions, the ones with the most realistic density 
and temperature profiles are those with rig « 0.05 cm~'^ 
within 10 kpc, because the density profiles of solutions 

^ Astute readers may recognize that this expression for the 
Field length as a function of entropy is somewhat larger than in 
Donahue et al. (2005). This d iscrepanc y arises from the normal- 
ization of the cooling function. IDonahu"e et al.i (i2005i) used a value 
of A appropriate for solar-metallicity plasma at ~ 1 keV, which is 
dominated by line cooling and cannot be extrapolated to higher 
temperatures according to A oc T^'^. This paper normalizes the 
cooling function so that it is appropriate for higher-temperature 
plasma, meaning that the assumed T^'-^ dependence leads to an 
underestimate of A, and therefore an overestimate of Ap, for tem- 
peratures ~ 1 keV. 
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Fig. 3. — Conductively balanced solutions assuming fc = 1/3 
for the intracluster medium of a 5 X lO^** Mq cluster. Each solu- 
tion depends on the temperature and electron density at the inner 
boundary, which are T/Ta = 0.5 (solid lines) and 0.3 (dotted lines) 
and Tie = 0.45, 0.15, and 0.045 in cm"''. Solutions with greater 
rte profiles (upper panel) have hotter temperature profiles (middle 
panel) because a greater temperature gradient is needed to off- 
set cooling. However, they all have similar entropy profiles (lower 
panel), which are close to the analytical approximation given by 
equation II20I I and shown by the dashed line. 



with larger inner densities are too shallow and represent 
a total amount of gas that is far too large. However, 
despite their differences, all of these solutions trace ap- 
proximately the same locus in the K-r plane, which is 
quite similar to the analytical approximation found by 
combining equations (fT6|) and (fT9l) with ar = 0.3 and 
a. = 5/2: 
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-1/3 
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This entropy profile, illustrated with a dashed line in the 
bottom panel of Figure [31 is therefore the characteristic 
entropy profile one expects of cool-core clusters in which 
electron thermal conduction balances free-free cooling. 
Its power-law slope in radius is a = 2/3, and its normal- 
ization depends on the value of fc, independent of cluster 
temperature, as long as A(r) oc T^/^. Above this line, 
conduction can transfer heat more quickly than radiative 
cooling can shed it. Below this line, conduction cannot 
compete with cooling. 

In groups and low-mass clusters, the characteristic en- 
tropy profile for conductive balance has a slightly greater 
normalization, because collisionally excited emission- 
line cooling increases A(T) and therefore the value of 

kTue at which conduct ion can balance cooling. Us - 
ing the cooling functions of lSutherland fc Dopital (|1993[ ) , 
we have derived a metallicity-dependent correction term 
that accounts for this additional cooling in systems with 
temperatures > 1 keV: 



K{r) « 5.0 rj 



2/3 / /( 



kpc \^-^/3 
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(21) 



The Z/Zq term in this expression is the gas metallicity 
in solar units. 

Figure |4] illustrates the approximations that equation 
1211 gives for conductively balanced profiles and compares 
them with the K{r) profiles of actual clusters. The thick 
solid line applies to clusters hotter than 3 keV, for which 
emission-line cooling is a small correction. The thick 
dashed line shows the locus of conductively balanced pro- 
files at lower temperature {^ 1 keV) and solar metallicity. 
Both lines are for /c = 1/3, and a scale bar shows how the 
normalizations of the lines would change for fc = 0.1 or 
1. Also shown are th e observed entropy profiles (from 
ICavagnolo et al.' '20 09D of two gro ups of clusters from 
the study of Raffcr tv et al.l (|2008[ ). Clusters with star- 
forming central galaxies (blue dashed lines) have K{r) 
profiles that generally descend below the solid line, into 
the regime in which cooling operates more quickly than 
conduction. Clusters with no star formation and no Ha 
in the central galaxy (red dotted lines) remain within the 
regime in which thermal conduction is more efficient than 
radiative cooling, at least for fc > 0.1. (The subset of 
clusters with Ha but without a color gradient indicating 
star formation in the central galaxy is not shown.) 

As discussed bv lVoit et al.l ( 2008() . these findings sug- 
gest that the critical locus for conductive balance divides 
galaxy clusters into two different types. Clusters whose 
entropy profiles remain entirely above this locus cannot 
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Fig. 4. — Comparison of observed intracluster entropy profiles 
with the loci of conductively balanced profiles. Thick lines show 
analytical approximations from equation II21I I for conductively bal- 
anced clusters with temperatures greater than 3 keV (solid line) 
and for conductively balanced groups with feT Ri 1 keV and solar 
metallicity (long-dashed line). Both lines assume /c = 1/3, and the 
scale bar shows how the vertical positions of both lines would shift 
for /c = 1/10 and 1. Thin lines show observ ed entropy profiles 
for the clusters discussed in IVoit et aTl II2008I1 . Profiles of clus- 
ters with star-forming central galaxies are illustrated with (blue) 
dashed lines. Those of clusters without star formation or Ha emis- 
sion in the central galaxy are illustrated with (red) dotted lines 
and remain in the regime in which thermal conduction can offset 
radiative cooling. 



maintain a multiphase intracluster medium because con- 
duction will evaporate clumps of gas that are cooler and 
denser than the ambient medium. With lower gas densi- 
ties and longer cooling times, the cores of these clusters 
are also more susceptible to dynamical heating, turbu- 
lent heat transport, and mergers. On the other hand, 
multiphase gas will tend to linger in the cores of clusters 
whose entropy profiles drop below the conductive locus, 
allowing them to harbor emiss ion-line nebulae and on - 
going star formation (see also iNipoti fc Binneyl 12004 ) . 
Because of their low central entropy levels, these clusters 
have cool, dense cores that are more difficult to disrupt. 

Conductive balance in this context should not be con- 
fused with thermal stability. Conductively balanced solu- 
tions for cluster cor es are well known t o be unstable (e.g . 
iBregman fc DavidI [1988; Soker 200^ iGuo et all l200a ). 
Conversely, environments that are out of conductive bal- 
ance are not necessarily prone to local thermal instabil- 
ities and growth of multiphase structure. This means 
that the critical locus in Figure |4] is more meaningful as 
an indicator of the need for thermal feedback, in addition 
to the heat supplied by conduction, than as a threshold 
for local thermal instability. 

A closer look at the entropy profiles in Figure |4] sup- 
ports the idea that the locus of conductive balance is a 
critical threshold for triggering thermal feedback. Above 
that locus, the entropy profiles shown by the dashed lines 
have a ^ 1 — 1.2 slopes expected of both gravitational 
heating and pure cooling. But below that locus these en- 



tropy profiles become shallower, with slopes that resem- 
ble the a — 2/3 slope of the locus itself. One way to un- 
derstand this behavior is to consider what happens to a 
cluster when conduction can no longer balance core cool- 
ing. As cooling proceeds, entropy levels in the core will 
drop and the inward mass flux will rapidly increase from 
zero toward the prodigious pure-cooling rate of equation 
(fT4|) . If thermal feedback does indeed depend on ac- 
cretion of intracluster gas at the cluster's center, then 
it should achieve a quasi-steady balance with radiative 
cooling with an entropy profile not far below the locus of 
conductive balance. This may account for the similarity 
in core entropy, density, and temperature profiles among 
cool-core clusters. 

5. THERMAL FEEDBACK 

Most of the recent attempts to explain how cluster 
cores stave off catastrophic cooling have focused on feed- 
back from an active galactic nucleus in the central galaxy. 
The heating per unit volume caused by feedback can be 
expressed with a term H. that depends on the mass ac- 
cretion rate at the center of the cluster. If conduction 
is also operating, the steady-state entropy equation then 
becomes 

P V dXviK ^, o. Id, dT .. 

7--T^ — ^^-p^^-y^^-r 22 

7 — 1 r d In r r'^ dr dr 

The right hand side of this equation must be negative 
everywhere in order to obtain a solution with steady in- 
flow (u < 0) and a positive entropy gradient. Even if 
the accreting matter condenses into a cold phase before 
feeding the active nucleus, equation ((22|) will remain valid 
outside the region where the condensation rate is a sig- 
nificant fraction of the mass inflow rate. 

This section uses equation ([^ to explore steady-state 
core conflgurations in which conduction, radiative cool- 
ing, and thermal feedback are all significant. It flrst 
presents a simple formalism for expressing the radial dis- 
tribution of thermal feedback and its dependence on the 
mass accretion rate. It then attempts to infer the ra- 
dial distribution of thermal feedback from observations 
of actual clusters, under the assumption that they are 
in a steady state with fc = 1/3. The main findings of 
this section are that thermal conduction may often be 
the most important source of heat in a cluster's core and 
that thermal feedback, when it is necessary, is rarely re- 
quired beyond ~ 30 kpc from the cluster's center. 

5.1. A Simplistic Model for AG N Feedback 

Thermal feedback from a central AGN can be phe- 
nomenologically modeled by expressing the outward en- 
ergy flux available for heating as 
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where e is an efficiency factor for thermal feedback and 
h{r) is a function describing how that energy flux declines 
with radius as it is converted into heat energy within the 
intracluster gas. The corresponding volumetric heating 
rate is 

„ , , , eMc^ dh ,„ , , 

^ ' 47rr2 dr ^ ' 

According to the sign convention adopted here, both M 
and dhjdr are negative when AGN feedback is operating. 



Grouping the terms proportional to M facilitates a 
direct comparison of heat advection with thermal feed- 
back: 
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The quantity in square brackets is always positive for 
clusters with a rising entropy gradient, and the relative 
importance of the two terms within the brackets de- 
pends on how the thermal energy of the gas compares 
with its rest-mass energy times the efficiency factor e. 
For a typical ICM temperature of ~ 5 keV, the magni- 
tude of the advection term within the square brackets is 
~ kT/finipC^ ~ 10~^. A value of e significantly exceeding 
10^^ is therefore necessary for feedback to suppress cool- 
ing and infiow. When this is the case, the mass-inflow 
rate needed to maintain a steady state is inversely pro- 
portional to e, and as long as the configuration is close 
to hydrostatic equilibrium, steady-state solutions with 
identical values of eM will be indistinguishable. Fur- 
thermore, the presence of a non-negligible heating term 
in the steady-state solution allows the entropy at radii 
where heating is significant to be smaller than in a con- 
ductively balanced solution. 

5.2. Inferring the Radial Distribution of Thermal 
Feedback 

The responsiveness of steady-state entropy profiles to 
the radial distribution of heat input raises the possibility 
that one might be able to infer the radial distribution 
of feedback heating in real clusters from observations of 
their temperature and density profiles. If a cluster's core 
configuration is truly steady, the observed density and 
temperature profiles can be plugged into the following 
equation to determine the steady-state distribution of 
heat input: 



dh 
dr 
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2 , 1 d 2 dT 
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Unfortunately, cluster temperature profiles are often too 
noisy for differentiating twice, so instead we will take 
advantage of the hydrostatic model of § [2] It gives tem- 
perature profiles in terms of a gravitational potential 
represented by T^(r) and an entropy profile shape. In 
what follows, the potential is assumed to be a simple 
NFW potential parametrized by Ta and ca, and the en- 
tropy profile is parametrized using the fitting formula of 
[Cavagnolo et al. (2009): K{r) = i^o + /'i:ioo(?'/100kpc)". 
For most of the clusters considered here, this five- 
parameter hydrostatic model adequately describes the 
observed temperature and density profiles. 

Figure [5] shows an example of a cluster that may be 
close to conductive balance. The bottom three panels 
show the model used to infer the radial distribution of 
thermal feedback in Abell 1795 and how it compares with 
the observed density, temperature, and entropy profiles 
of this cluster from the compilation of ICavagnolo et aLl 
()2009( l. The top panel shows the relative importance of 
heating, cooling, and conduction as a function of radius, 
assuming that this cluster is in a quasisteady state with 
negligible advection of heat. Under this assumption, the 
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Fig. 5. — A quasi-steady solution for Abell 1795. The top panel 
shows the relative importance of cooling (solid line), conduction 
(long-dashed lines), and thermal feedback (dotted line) in terms of 
timescales. The lower a line is in this panel, the more rapid and 
important the process. Solid lines in the lower three panels show 
the density, temperature, and entropy profiles, respectively, of the 
solution used to compute the cooling and conduction rates from 
which heating was inferred. Dashed lines in those panels show 
the observations and shading shows uncertainty ranges for these 
quantities from the database of lCavagnolo et al.i 1,2009,) . 
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heating rate can be inferred from equation (j22l) by set- 
ting the left-hand side to zero and plugging in the model 
profiles of temperature and density. A suppression fac- 
tor fc = 1/3 has been assumed. A solid line in the top 
panel shows the cooling time tc, defined as the quotient 
of the heat content at constant pressure, 5pkT/2fj,mp, 
and the radiative cooling rate. The dashed lines (blue 
and green) show the conduction time icond, equal to the 
quotient of the heat content and the conduction term in 
equation (P^ . The dotted line (red) shows the inferred 
heating time th determined by dividing the heat content 
by the difference of the cooling and conduction terms. 

Within 100 kpc, the timescales for cooling and con- 
duction in Abell 1795 are nearly equal, meaning that the 
conduction and cooling terms in e quation (|22l) nearly 
offset each other (see also iZakamska fc Naravan 2003) . 
Thermal feedback is therefore not necessary to keep this 
system in a steady state, and the inferred heating term 
corresponds to a long timescale (> 10 Gyr). Outside 
of 100 kpc, where the temperature profile peaks, con- 
duction acts as a coolant channeling heat away from the 
peak. The line representing icond is green in this region 
instead of blue to indicate the change in sign of the con- 
duction term. Under the steady-state assumption, heat- 
ing must compensate for conductive cooling here, and the 
red dotted line for th indicates the implied heating rate. 
However, the timescales in this region are comparable to 
the age of the universe, so the steady-state assumption 
is unlikely to be valid. In fact, mergers and other forms 
of dynamical heating not included in equation (|22p are 
relevant on these timescales. We will therefore restrict 
our attention to radii < 100 kpc, where the steady-state 
assumption has a greater chance of being applicable. 

Figure |6] shows a different example, the cluster 
Abell 1835. The lines represent the same things as in 
Figure [5l but the implications are different. Thermal 
conduction in this system cannot offset radiative cool- 
ing at small radii. The inferred thermal feedback rate 
within that region is greater than the conductive heating 
rate, causing the heating timescale t^ shown in the figure 
(dotted line) to be shorter than the conductive timescale 
icond (dashed line), but only within < 30 kpc of the clus- 
ter center. 

In that respect, Abell 1835 is representative of most 
of the clu sters with star-formi ng central galaxies in the 
sample of Raf fertv et al.l ()2008f) . Figure [7] shows the ratio 
of inferred heating to cooling as a function of radius in all 
18 of those clusters. Only three of those systems (dotted 
lines) are close to conductive balance at all radii. An- 
other 11 (solid lines) require thermal feedback in order 
to remain in a steady state, at least for fc = 1/3, but in 
most of those cases inferred feedback heating dominates 
over conductive heating only within < 30 kpc. There are, 
however, four systems in which the inferred steady-state 
heating rate remains large out to 100 kpc. Part of the 
reason for the rising level of implied heating near 100 kpc 
in these systems is that the temperature profile is reach- 
ing its maximum value there, causing a decline in the 
conduction rate inferred from the temperature gradient. 
It is possible that other forms of diffusive heat transport, 
such as turbulent transport, could be more important 
than conduction in these regions, reducing the need for 
thermal feedback there. 

Integrated over the whole cluster core, conduction may 
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Fig. 6. — A quasi-steady solution for Abell 1835. The top panel 
shows the relative importance of cooling (solid line), conduction 
(long-dashed lines), and thermal feedback (dotted line) in terms of 
timescales. The lower a line is in this panel, the more rapid and 
important the process. Solid lines in the lower three panels show 
the density, temperature, and entropy profiles, respectively, of the 
solution used to compute the cooling and conduction rates from 
which heating was inferred. Dashed lines in those panels show 
the observations and shading shows uncertainty ranges for these 
quantities from the database of lCavagnolo et al.i 1,2009,) . 
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Fig. 7. — Inferred ratio of thermal feedback to radiative cooling 
as a function of radius for conductive clusters (fc = 1/3) assumed 
to be in a stead y state. Thes e clusters are the 18 cool-core clus- 
ters from Raffcr tv et al.l II2008I ) with star-forming central galaxies. 
Dotted (blue) lines show the three clusters that are nearly in con- 
ductive balance. Solid lines show the 11 clusters in which con- 
duction is more important than thermal feedback outside of the 
central ~ 40 kpc. Dashed (red) lines show the four clusters in 
which thermal feedback is needed to maintain a steady-state in the 
outer parts of the cluster core. 



be a more important heat source than thermal feedback. 
Figure |8] compares the volume integral of the imphed 
heating rate within a radius of 100 kpc to the integrated 
coohng rate. In most cases the inferred amount of ther- 
mal feedback is less than 50% of the total cooling, mean- 
ing that conduction with fc = 1/3 can replace the ma- 
jority of the heat losses in these systems. There does, 
however, appear to be a trend in the upper envelope of 
this relationship. At lower temperatures there are more 
objects in which thermal feedback seems to be the dom- 
inant heat source. It would appear that conduction is 
always important in the highest-temperature cool-core 
systems but is less effective at compensating for cooling 
in at least some of the lower-temperature cool-core sys- 
tems. 

Taken as a whole, these result s are in alignment with 
the fin dings oflVoigt et al.l (I2002D. Zakamska fc NaravanI 
(|2003D . and |WgL&Fabian| ^QM). Those studies also 
found that thermal conduction could plausibly balance 
radiative cooling in some cluster cores, but that it fails to 
balance cooling within the central ~ 30 kpc in many oth- 
ers, necessitating a feedback response to prevent cooling 
rates in excess of the observed limits. This study differs 
from those earlier efforts in that it fits a physically mo- 
tivated hydrostatic model to the cluster data in order to 
obtain a numerically differentiable temperature profile, 
which it then uses to infer the radial distribution of ther- 
mal feedback. Those distributions of heat energy are un- 
likely to be accurate in detail, because cool-core clusters 
are probably not in a perfectly steady state, and because 
isotropic conduction with fc — 1/3 is probably overly 
simplistic, but they do indicate that AGN feedback does 
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Fig. 8. — Inferred ratio of thermal feedback to radiative cool- 
ing integrated within 100 kpc for the same clusters as in Figure [7] 
Squares show clusters close to conductive balance. Circles show 
clusters in which conduction is more important than thermal feed- 
back outside of the central ~ 40 kpc. Triangles show the four 
clusters in which thermal feedback is needed to maintain a steady- 
state in the outer parts of the cluster core. 

not necessarily need to match radiative cooling all the 
way out to ~ 100 kpc. Heat input in the vicinity of the 
observed X-ray cavities may be all that is necessary. 



6. SUMMARY 

Inspired by the fact that the radial profiles of intra- 
cluster gas density, temperature and entropy in cool-core 
clusters bear a strong resemblance to one another, this 
paper has tried to identify the physical processes that 
govern the form of those profiles. Because the similar- 
ity among these profiles extends to within ~ 10 kpc of 
the cluster's center, where cooling, heating, and dynam- 
ical processes can operate on timescales < 100 Myr, we 
have assumed that cool-core clusters have settled into 
a quasi-steady state that is close to hydrostatic equilib- 
rium. Proceeding from that assumption, the paper de- 
rived the general characteristics of various quasi-steady 
configurations for the intracluster medium. 

First, it showed in § 2 that the temperature profiles of 
cool-core clusters simply reflect the shape of the under- 
lying gravitational potential, as long as the intracluster 
entropy profile has a moderately positive radial gradient. 
When that is the case, the gas temperature at each radius 
is determined primarily by the value of T^ oc M{< r)/r 
at slightly larger radii and to a lesser degree by the slope 
of the entropy profiles. The temperature does noi depend 
on the absolute value of gas entropy — only the slope of 
the entropy profile matters. For an NFW-like potential 
with a concentration ca ~ 4, the gas temperature peaks 
in the O.Ita — 0.2rA range, and within a given potential 
well, configurations with shallower entropy slopes peak at 
greater temperatures and smaller radii. Adding a small 
isothermal potential at the center to represent the mass 
of a central cluster galaxy helps to explain why clusters 
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contain so little gas at temperatures less than 30% of the 
virial temperature. 

Next, § 3 derived a set of pure-cooling solutions in 
which inward advection of heat balances radiative cool- 
ing. Steady cooling of gas in a cluster potential turns 
out to produce entropy profiles similar in slope to the 
K{r) oc r^-^ profile generated by gravitational structure 
formation. This coincidence accounts for why there is no 
break in slope in the entropy profiles of cool-core clusters 
at radii where the cooling time is similar to the age of 
the universe. Furthermore, the mass-inflow rates of these 
steady pure-cooling solutions are oc K~^ and greatly ex- 
ceed the observed star-formation rates of central galaxies 
in cool-core clusters. This implies that entropy profiles 
in cool-core clusters should never drop far below the crit- 
ical profile at which inflow begins to fuel feedback. Once 
inflow begins, the feedback response should rapidly rise 
as the entropy level at a given radius drops. 

Section 4 derived the characteristics of solutions that 
are in conductive balance, showing that despite widely 
differing density and temperature proflles, steady conflg- 
urations with the same value of fc trace the same locus in 
the K(r)-r plane. Simple expressions for this locus are 
given for hot systems (> 3 keV) and for cooler systems 
in which the cooling rate is metallicity dependent. Com- 
paring the observed entropy profiles of galaxy clusters to 
these loci shows that clusters without star-forming cen- 
tral galaxies or other evidence of multiphase gas all lie 
above the locus for conductive balance, suggesting that 
conduction is preventing radiative cooling from produc- 
ing a multiphase intracluster me dium that can a ccrete 
into the central galaxy (see also IVoit et al.ll2008l ). On 
the other hand, the entropy profiles of classic cool-core 
clusters tend to drop below this critical locus at small 
radii, but they do not drop very far below it. Instead, 
their entropy slopes bend in the vicinity of the critical lo- 
cus and tend to track it at slightly lower entropy levels, 
suggesting that AGN feedback is triggered when an en- 
tropy profile drops below the critical locus and prevents 
the entropy from dropping much further. This would ex- 



plain why the density, temperature, and entropy profiles 
of cool-core clusters are so similar to one another. 

Section 5 considered cluster configurations in which 
cooling, conduction, and thermal feedback all have roles 
to play in maintaining a steady state. Instead of assum- 
ing a particular model for thermal feedback, the paper 
attempted to infer the spatial distribution of feedback 
heating from observations of cool-core clusters. These 
were fit to a five-parameter model in which the gravi- 
tational potential has an NFW form and the intraclus- 
ter entropy profile is a power law plus a constant. This 
model provides a smooth temperature profile that can 
be differentiated to obtain the heat deposited by ther- 
mal conduction as a function of radius. If heat advection 
is assumed to be negligible, the steady-state heating rate 
is then equal to the difference between radiative cool- 
ing and conductive heating. Some conductive cool-core 
clusters do not require thermal feedback to offset cool- 
ing, but most of them do. However, the inferred thermal 
feedback rates usually exceed conductive heating only 
within < 30 kpc, where X-ray cavities are most often 
found. Conduction with /c = 1/3 can account for most 
of the heating required beyond that radius, and when in- 
tegrated over the entire core out to 100 kpc, conduction 
is a more important heat source than thermal feedback 
in most of the clusters considered here. If conduction 
is indeed as important as these results imply, then it is 
essential to include it in numerical simulations that seek 
to accurately model the core structures of galaxy clus- 
ters, including the stellar content of the brightest cluster 
galaxy 
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